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RESIDUES  OF  INTEGRALS  WITH  THREE-DIMENSIONAL  MDLTIPOLE 
SINGULARITIES ,  WITH  APPLICATION  TO  THE  LAGALLY  THEOREM 

I,  INTRODUCTION 

An  important  mathematical  relation  in  the  theory  of  three-dimen¬ 
sional  irrotational  flow  is  the  Gauss-Green  transformation  between  volume 
and  surface  integrals,  A  very  useful  result  of  this  transformation, 
Green's  third  formula,  requires,  in  its  derivation,  the  evaluation  of  the 
limit  of  a  singular  integral  over  the  surface  of  a  sphere  as  the  radius 
of  the  sphere  approaches  zero.  Since,  in  this  case,  the  singularity  is 
due  to  the  potential  of  a  source  at  the  center  of  the  sphere,  its  limit 
may  be  called  the  residue  of  a  source.  Similarly,  limits  of  integrals 
over  the  surface  of  a  sphere  of  vanishingly  small  radius,  with  higher- 
order  derivatives  of  the  source  potential  in  the  integrand,  will  be 
called  residues  of  multipoles •  The  latter  occur  in  the  derivation  of  the 
Lagally  theorem  for  the  force  and  moment  acting  on  a  body  moving  in  an 
irrotational  flow  when  multipoles  are  present  in  the  hydrodynamic 
singularity  system  within  the  body  (see  Landweber,  1967  and  Landweber  and 
Miloh,  1980,  1981). 

In  contrast  to  the  very  simple  derivation  of  the  residue  occurring 
in  Green's  third  formula,  the  evaluation  of  the  multipole  residues  was  a 
challenging  application  of  the  theory  of  spherical  harmonics.  The 
derivations  of  a  set  of  multipole  residues,  which  were  required  but  not 
included  in  the  aforementioned  references,  will  be  presented  here. 

This  work  is  dedicated  to  my  long-time  friend  and  colleague,  Chia- 
Shun  Yih,  with  whom  I  worked  30  years  ago  on  the  generalization  of  the 
Lagally  theorem  to  unsteady  flows  (1956).  Support  by  the  Office  of  Naval 
Research  under  Contract  N0001 4-83-K-01 36  (NR-062-183)  is  gratefully 
acknowledged. 


II.  FORMULATION  AND  TABULATION  OF  RESULTS 

The  velocity  potential  $  of  a  multipole  of  strength  M  and  of  order 
<1  -  a+8+Y/  si-tuated  at  a  point  p(xs,ys,zs),  with  position  vector  r  ,  of 
a  rectangular  Cartesian  coordinate  system  (x,y,z)  will  be  expressed  as 


1 


*P  =  -  mds(R>'  Ds  = 


9x  9y  9 
s  s  s 


— ,  R2  =  (x-xs)2+(y-ys)2+(z-Zs)2 


(1  ) 


Let  S  denote  the  surface  of  a  sphere  of  radius  c  with  center  at  the  point 


-  r 


It  is 


0(x  ,  y  ,  zn)  and  position  vector  r  such  that  c  >  b  =  r 

o  o  O  q  II 

necessary  to  displace  the  singularity  from  the  origin  in  order  to  ensure 
that  the  coordinates  (x,y,z)  of  a  point  on  S  are  independent  of  (xs,ys/ 

Zg ) •  PUt 


p  =  r  -  r  =pn,  p  =  r  -  r  =bn  (2) 

o  p  p  o  p 

where  n  is  the  unit  vector  in  the  direction  from  0  to  P  and  n  the  unit 
P 

vector  from  0  to  the  point  (x,y,z). 

The  integrals  of  interest  involve  a  pair  of  vectors  u^  and  v  where 

u  =  V$  t  v  =  V¥,  and  ¥  is  regular  harmonic  within  the  sphere  and  on 

P  P  __  __ 

S.  Thus  both  u  and  v  are  irrotational  and  solenoidal  within  the  sphere 
P 

and  on  S,  except  for  the  multipole  singularity  of  at  P.  Results  for 

the  residues  are  given  in  Table  1 •  Their  derivations  will  be  presented 

in  the  following  sections. 

We  observe  that,  when  q  =  0,  the  first  integral  in  Table  1  reduces 
to  the  residue  of  Green's  third  formula, 

lim  f  ¥  ri«V(1/R)dS  =  -  4ttM^(x  ,y  ,z  )  (3) 

X  s  s  s 

c+o  S 

The  second  integral  is  a  corollary  of  the  first  since  each  of  its 

components  is  of  the  type  of  the  first.  This  shows  that  the  second 

integral  does  not  require  the  conditions  Vxv  =  0  and  V*v  =  0.  it  is 

included  in  the  list  of  integrals  because  pairs  of  the  three  integrals, 

(5),  (6)  and  (7),  yield  the  residues  for  the  three  vector-product 

combinations  of  u  ,  v  and  n,  viz: 

P 
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Table  1 .  Multipole  Residues 


lira  f  n-Up^dS 

C  -HD  S 


4  ttM  ( q+ 1 ) 

2q+1 


q  * 
°>>P 


(4) 


/  "-“p  “ "  -  i^sl  »;«*», 
c+o  S 


(5) 


lim  /  Up  *v  ndS 
040  S 


_ 4  ttM _ 

(2q+1 ) (2q+3) 


Dj(v)p 


(6) 


lim  /  v«n  u  dS  =  ^ ^  D^(v) 

1  P  2q+3  s  P 

C-KD  S 


(7) 


lim  /  n.uppxvdS  =  [D^(rxv)p  -  rpxD^(v)p]  (8) 

C  -K)  S 


lim  f  v.n  pxupdS  =  [D^(rxv)p  -  rpxD^(v)p]  (9) 

C  40  S 


lim  f  $pv-ndS  =  "|^Dg(T)p  <10) 

C40  S 


lira  f 
c+o  S 


$  ndS 
P 


M<5 


(11) 


lim  f  $prxndS  =  -  M{  rpxD^  rp  (12) 
c-vo  S 


*  /  \ 

(  )p  indicates  that  the  quantity  between  the  parentheses  is  evaluated  at 
the  point  P(xs,ys,zs). 
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lim  /  nx(vxup)dS  =  T2q+!)|lq+3)  Ds(v)P  (13) 

C-HD  S 

lim  /  vx(upxn)dS  =  D^(v)p  <14> 

C-K)  s 

lim  f  upx(nxv)  dS  =  "  ^q+f^  Ds(v)P  (15) 

C  -*0  S 

III*  SOME  LEMMAS 

It  will  facilitate  the  derivations  if  we  first  collect  some  useful 
spherical-harmonic  properties.  Put  £  =  x-xq,  r\  =  y-y  ,  £  =  z-zQ/  so 
that  p  has  the  components  (£,n,£).  We  shall  also  employ  spherical 
coordinates  w  ~  cos  9,  with  the  origin  at  the  point  0  and  the 

polar  axis  along  OP. 

a)  The  scalar  Y  and  the  vector  v(£,g,£)  may  be  expanded  in  the 

forms 


¥  =  E  ¥  UfTlfCJ*  ¥  =  pmS  (p,<)>)  (16a) 

^  m  m  m  T 

m=0 


v 


E  v  ( £,n, O ,  v 
~  m  m 


p"1  T  (  y ,  cj) ) 
m 


(16b) 


where  ¥  and  v  are  solid  spherical  harmonics,  homogeneous  of  degree  m 
m  m  _ 

in  5,g,£,  and  Sm  and  T  are  the  corresponding  surface  spherical  harmonics 
of  order  m. 

b)  Since  ¥  and  v  are  homogeneous  in  £,ri,£,  they  satisfy  the  Euler 
in  m 

relations 


p*V¥  =  m'F  ,  p»Vv  =  mv 
mm  mm 


(17) 
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c)  When  the  irrotational  vector  v  is  regular  harmonic,  the  quanti¬ 
ties  pxv  and  p*v  are  also  regular  harmonic. 


Proof ;  In  index  notation  with  the  summation  convention,  we  have 


2 

3  v. 


3v% 


35  Hjk  k 


5-ivv  “  eijk  9^ 

£  1 


+  2 


3E  6j  9}  2eijk  35. 


=  0 


3v. 


-  E.  v. 

3£  3E  H  i 


3  v.  3v. 

_ 1 1  X 

3?„35„  +  2  liTT  =  2  ~h7 


=  0 


where  e  is  the  permutation  tensor,  and  6  is  the  kronecker  delta 
iDk  ij 

function. 

d)  The  orthogonality  relation  between  Legendre  polynomials  P  (y) 

_  n 

and  a  surface  spherical  harmonic  Sm  or  T  shows  that  the  integral  of 

their  product  over  the  unit  sphere  is  proportional  to  the  value  of  Sm 

or  T  when  y  =  1;  (see  MacRobert,  1945).  Since  Sm  or  T  is  homoqeneous 
m  m  m 

of  degree  m  in  £,n,£,  and  the  point  P  lies  on  the  polar  axis,  this  result 

may  be  expressed  in  terms  of  ( ip  )  =  p  (S  )  or  (v  )  =  pm(T  )  , 

mPPmP  mPpmP 


e) 


1  2  rr  S 

bm  /  /  Pn(u)l_m}dpd<|, 

-1  o  T 

m 


4tt6 

mn 

2n+1 


(¥  ) 

{  _m  P) 
(v  ) 
m  P 


(18) 


1  2tt 

h11"1  /  /  Pn(  lj)e"L^Sm(  U'  <t))dUd<t> 

-1  0 


4tt6  3¥  S'? 

mn  _ m  .  _ m> 

2n+1  (  3y  +1  3z  ;p 


(19) 


where  P’(y)  is  the  associated  Legendre  function  of  degree  n  and  order  1, 
_ n 

i  =  /-I,  and  the  (x, y, z ) -coordinate  system  has  been  taken  so  that  the  x- 
axis  is  parallel  to  the  polar  axis  of  the  spherical  coordinate  system. 

Proof:  S  (  y ,  c{> )  may  be  expanded  in  the  series 

m 

S  =  AP  (y)  +  P^(y)  (B  cos  d>  +  C  sin  d> )  +  ... 
m  m  n  T  T 
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Then,  by  orthogonality,  we  have 


I  f''l  <>.)eiVl.d*  ‘  2gffi±L  Vn<B+ic)  <201 

-1  0 


In  terms  of  unit  vectors  e^  ,  e^,  e^  in  the  directions  of  increasing 
values  of  (p,0,((>),  the  gradient  of  Y  at  (1,0,<j>)  is 

(  K  ^  x  =  ei A  +  T  m(m+1  )  [e0  (B  cosd>  +  C  sin  6)  +  e0  (-B  sin(J>  +  C  cosd)  ] 
ml,0,c})l2  z  T  T  z  T  T 

Then,  taking  components  of  the  gradient  in  the  y-  and  z-directions ,  we 
obtain 

3y  1,0, (j)  =  2  (9z  ^  1  ,0,4,  =  2  m(m+1)c 

which,  substituted  for  B  and  C  in  (20),  yields  (19),  the  result  we  wished 
to  prove. 

y 

rr  \li  /r  tti 

(21) 


lira  Dq  {1}  =  Dq  (_m}  5 

s  —  s  —  qm 
c  *0  v  v  ^ 


m 

since  the  terms  of  the  expansion  (16a,b)  with  m<q  vanish  under  q 
differentiations,  and  those  with  m>q  vanish  as  c-K). 

g)  Two  different  expansions  of  the  velocity  field  of  a  multipole 

u  will  be  useful.  We  have 
P 


U  =  MDq(R/R3),  R  =  p-p 
p  s  P 


(22) 


The  expansion  of  1 /R  then  gives 


R  1  00  bn  _  ,  x  - — 

—  -  R  -  7s  E  -JFT  p„(>,)-  »  -  "-n"'  n- 

R  n=1  p 


P'  P 


(23a) 


and 


“  n“1 

=  l  — —  [p  (|i)n  -  ]?  (y)n  ] 
3  u  n+1  n  n-1  P 

R  n=1  p 


(23b) 


6 


since 


V  (bnP  )  =  nbn_1P  n  +  bnP  V  (r^)  .  n 

s  n  n  p  n  s  b 


=  bn  1  [nP  n  +  £  (n-yn  )  ] 
n  P  n  P 


and  (see  MacRobert,  1947) 


yPn(y)  -  nPn(y)  =  Pn_-,(y) 


(24) 


Here  the  dot  indicates  differentiation  with  respect  to  the  argument.  We 
shall  also  require  the  relation 

Pn(y)  -  VPn_.,(u)  =  nPn_1(U)  (25) 

IV.  DERIVATIONS  OF  RESIDUES 

The  residues  in  Table  1  will  be  derived  in  the  order  shown,  except 
that  we  must  first  verify  (15)  in  order  to  derive  (6). 

A*  Derivation  of  (4)  and  (5).  Substituting  (16a)  for  y,  (22)  and 

(23b)  for  u  ,  putting 
P 

2 

dS  =  -c  dydcj) 


and  applying  (25),  we  obtain 


{  n*u  f  dS  =  MDq  Z  l 
S  P  S  m=o  n=1 


1  2  7T 

f  /  cm“n+1bn"1nP  (y)S  dyd (J> 
„  ^  n-1  m 


Orthogonality  of  the  surface  spherical  harmonics  requires  n=m+1 ,  so  that 
the  right  number  of  this  equation  becomes 


7 


or,  by  (18), 


MDq  l  (m+1)b 
s 

m=o 


1 


I 


2  7T 

/  P  (y)S  dyd<j> 
A  m  m  T 


lim  M 
c-kd 


47r(m+1  ) 

Z  2m+1 
m=o 


D ?(»)_ 
s  m  P 


47TM(q+1  ) 
2q+1 


Dq(  'P) 
s 


p 


since  the  terms  of  the  expansion  of  V  in  (16a)  vanish  under  q  differen¬ 
tiations,  and  those  with  m  >  q  vanish  as  c-K).  This  is  the  result  given 
in  (4)*  As  has  already  been  remarked,  the  residue  (5)  is  an  immediate 
consequence  of  (4). 

B#  Derivation  of  (6)  and  (15).  Substituting  (16b)  for  v,  and  (22) 

and  (23a)  for  u  ,  we  obtain 
P 


\  u  x(nxv)dS 
S  P 


MDq 

s 


E  t  c] 


.m-n+1 


1  2tt 


m=o  n=o 


bnV  x/  / 

S  -1  0 


P  (y)nxT  dyd((> 
n  m 


According  to  Lemma  c,  nxT  is  a  surface  spherical  harmonic  of  order  m+1 , 

m 

and  hence,  by  orthogonality,  the  integral  vanishes  unless  n=m+1 .  Then, 
by  (18), 


lim  /  UpX (nxv)dS  =  MD*VgX  ? 
c-K)  S  n=0 


4ttM 
~  2q+3 


Dq[  V  x  (  pxv  )  ] 
s  s  q 


P 


since  the  differential  operator  DqV  is  of  order  q  +  1 

s  s 


But 


V  x ( pxv ) 
s  P 


(v  *V  p  +  pV  «v 
q  s  s  q 


p  •  V  v 
s  q 


v  V 

q  s 


(q+2) 


since 


8 


and,  by  (17) 


v  *Vp  =  v  ,  V»v  =  0,  V*p  =  3 

q  q  q 


p  •  Vv  = 

q 


% 


Then,  by  (21 ) , 


lim  j  upx(nxv)dS 
c-kd  S 


4irM(q+2)  q,~  .  =  _  4TrM(q+2l  Dq,-) 

2q+3  s  q  P  2q+3  s  P 


as  we  wished  to  show. 

The  result  in  (6)  can  now  be  obtained  as  the  sum  of  the  residues  in 
(5)  and  (15),  since 


u  «v  n  =  u  x(nxv)  +  n  *u  v 
P  P  P 


C#  Derivation  of  (7).  By  the  same  procedure  as  in  the  previous 
cases,  the  integral  in  (7)  is  first  expressed  as  a  double  summation  and 
then  reduced  to  the  single  term 


lim 

c-kd 


4ttM 

2q+3 


DqV  (p*V  ) 
s  s  q  P 


since, 
=  VT 
by 


q+1 

17) 


by  Lemma  c 

where  'F 

q+1 


,  T  *n  is  a  surface  spherical  harmonic, 
m 

is  a  spherical  harmonic  of  degree  q+1 • 


Also 

Then 


put  v 


q 

we  have 


V  (  p  *  v  )  =  V  (p-VY  )  =  (q+1  )  V  (¥  )  =  (q+1  )  (v  ) 

s  q  P  s  q+1  P  s  q+1  p  q  P 


Hence,  by  (21 ) , 


lim  J  v«n  UpdS 
c-k>  S 


4”M(q+1)  nq(~  X  _  4 ttM( q+ 1)_  Dq(-} 
2q+3  sy  q;P  2q+3  s  P 


9 


as  was  stated 


D#  Derivation  of  (8).  Since  pxv  is  regular  harmonic  (Lemma  c) 
and  p=r-r  ,  we  obtain  from  (5) 


^ - -  .  4'irMq 

lim  J  n*Uppxv  dS  =  lim  2q+i 

c*o  S  C-K) 


Dq  [  (r-r  ) x v ] 
s  op 


4-TTMq  [Dq 
2q+1  s 


(r  x  v) 

P 


r  x  Dq  (v)  ] 
p  s  p 


in  agreement  with  (8). 

E#  Derivation  of  (9).  As  in  the  derivation  of  (7),  we  have 


p.v  =  p.V'P  =  £  p*Vf  =  E  m'P  =  E  mpmS 

mm  m 

m=o  m=o  m=o 


(26) 


Hence ,  by  (22)  and  (23b),  since  p=pn,  the  integral  in  (9)  becomes 


,n-1— 
oo  co  mb  n 


1  2tt 


I  =  /  v»n  pxu  dS  =  MDq  E  E  - Jp  x  f  f  n  P  .(yJS  dydrf) 

9  J  P  s  „  n-m-1  i  n-1  p  m  M  y 

S  m=o  n=1  c  -10 


We  have 


n  -  \xn  +  sin  0  (j  cos  cb  +  k  sin  6) 
P 


(27) 


where  n^,  j  and  k  are  an  orthogonal  set  of  unit  vectors.  Hence,  applying 

(19)  and  putting  bn  =  p  ,  we  obtain 

P  P 


pPx^vm-1^P  I9  "  MDs  S  2m+1  pPx(vm-1^P  2q+^ 

m 


Dq(  pxv) 

s  ^  p 


and  then,  since  p^  =  r^-r^, 


lim  1^ 

C  -HD 


[Dq  (rxv)  -  r  xDq  (v)  ] 

2q+1  s  P  P  s  P 
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in  agreement  with  (9)* 

F.  Derivation  of  (10).  The  desired  result  is  immediately  obtained 
by  writing 

n 

$  =  -  MDjf  E  “ P  (p)  (28) 

p  s  n+1  n  K 

n=o  c 

and  applying  (26) ,  (18),  and  (21). 

G.  Derivation  of  (11).  Since  n  is  also  a  surface  spherical 
harmonic,  application  of  (27)  and  (28)  yields,  by  orthogonality, 


$pndS  = 


4ttM 
3  U1 q^sq ' 


4ttM 


<5,  D  r __ 

1q  sq  P 


since  bn  =  p  =  r  -r  .  This  agrees  with  (11). 

P  P  P  0 

H.  Derivation  of  (12).  Substituting  r  =  r  +cn  into  (12)  and 

o 

applying  (11),  we  obtain 


lim  f  $  rxndS  =  lim  r  xf  $  ndS 

J  T>  r\  J  P 


C-KD  S 


C-K) 


4tt 


M5  r 
1  q  P 


xD^r* 
s  P 


in  agreement  with  (12). 


V.  EXAMPLE 


Force  on  a  Stationary  Sphere  in  a  Steady,  Slightly  Nonuniform  Flow 

A  sphere  of  radius  c  is  introduced  into  an  irrotational  flow  of  an 
inviscid  fluid.  We  shall  employ  a  rectangular,  Cartesian  coordinate 
system  (x,y,z),  and  a  spherical  coordinate  system  (R,y,<j>),  u  =  cos 9,  with 


polar  axis  along  the  x-axis  and  both  with  origin  at  the  center  of  the 
sphere,  so  that 

x  =  R  cos  0,  y  -  R  sin  0  cos  <J>,  z  =  R  sin  0  sin  <|>  (29) 


The  velocity  potential  of  the  undisturbed  flow  will  be  written  as 


2  1 

V  =  URy  +  R  S  (p,<|>),  S  =  A  P  (]i)  +  P  (y)(A  cost))  +  B  sin<f>) 
2  2  o  2  2  1  1 


+  P2(y)(A2cos  2(f>  +  B2  sin  2(f)) 


(30) 


This  corresponds  to  the  most  general  harmonic  flow  quadratic  in  x,y,z. 
The  x-axis  has  been  taken  in  the  direction  of  the.  uniform  stream, 
represented  by  the  first  term.  Its  expression  in  terms  of  x,y,z  is 

Y  =  Ux  +  ~  Aq(2x2  -  y2  -  z2)  +  3A^xy  +  3B^xz  +  3A2(y2-z2)  +  3B2yz 
2  12  12 

=  Ux+Aqx  +(3A2-“AQ)y  -(3A2+‘JAo)z  +  3B2yz+3B,j  zx+3A^  xy  (31) 

The  velocity  field  corresponding  to  y  is  then 

v  =  V¥  =  i ( U+2A  x+3A  y+3B  z)+j[3A  x+(6A  -A  )y+3B  z] 
o  111  2  o  2 

+  k  [3B  x  +  3B  y  -  (6A  +  A  )z]  (32) 

1  2  2  o 

where  i,j,k  are  the  unit  vectors  in  the  x,y,z  directions. 

The  disturbance  potential  $  which  satisfies  the  boundary  condition 


8R 


($  +  HO 


=  0 

R=c 


(33) 
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is 


3  5 

Uc  u  2c  _  ,  ,  v 

$  -  2  +  3  S2 <  Ur  <|>) 

2R  3R 


(34) 


Each  of  the  terms  of  (34)  can  be  expressed  in  terms  of  multipoles 
situated  at  the  center  of  the  sphere  by  applying  the  formula  given  by 
Hobson  (1931), 


-  — 


(n-m) ! 


n-m 

_3 _  ,_9 _  .  3  \  m  1 

n-m  1  3y  +  1  3z;  R 
3x 


(35) 


which  gives 


y_  =  _  _a_  j.  p2(„u.!.  =  iili 

2  3x  R'  3  2  2  R 

R  R  3x 


1  1 
P2(y)  cos  <})  g2  1  p2(u)  sin  <j»  g2  1 


R 


3x3y  R' 


R 


3x3z  R 


(36) 


P2(y)  cos  2 4>  g2  g2  1  ?2(  p)  sin  2$  g2  ^ 

3  =  (  2  "  2)  R'  3  =  2 

R  3y  3z  R 


3y3z  R 


The  force  on  the  sphere  can  now  be  obtained  from  the  Lagally  formula 
(Landweber,  1967) 


F  =  -  4irp  E  M  (D^v) 
s  s  o 
s 


Here,  by  (1)  and  (36),  the  successive  operators  are 

s 


(37) 
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2  2  2 

n1  =  JL_  2  t _ 9 _  _  _9 _ v 

s  3x'  Ds  “  ' 3x3y'  9x9z' * • * 9y 3z ; 


—  2 
Since  v,  given  in  (32),  is  linear  in  x,y,z,  the  second  derivatives  Dg  do 

not  contribute  to  F,  and  hence  F  becomes 


—  3  Jlv 

F  =  -  2 7T  pUc  (— — ) 

3x'o 


=  -  2TTpUc°(2iA  +  3jA  +  3kB 
o  1  1 


(38) 


We  observe  that  ~  Uc  is  the  doublet  strength  at  the  center  of  the 
sphere  in  a  uniform  stream.  The  result  in  (38)  is  the  well-known  Lagally 
expression  for  the  force  on  a  doublet  in  a  nonuniform  stream,  but, 
unexpectedly,  with  the  doublet  strength  unaffected  by  the  quadratic 
nonuni formi ty . 

The  foregoing  procedure  can  be  readily  extended  to  nonuniform  flows 
of  higher  order.  If  the  velocity  potential  contains  a  homogeneous 
harmonic  of  degree  n,  its  contribution  v^  ^  to  the  velocity  would  be  of 
degree  n-1 •  Then 


(Dq  V  )  =6  ,  (Dn  1  V  ) 

s  n-1  o  q,n-1  s  n-1  o 


(39) 


since  the  terms  with  q  <  n-1  vanish  at  the  origin.  The  force  (37)  would 
then  become 


F 


00 


4ttp 


I 

n=1 


M 

n 


(Dn  1  v  ) 
s  n-1  o 


(40) 
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Appendix 


An  Alternative  Procedure  for  Deriving  the  Multipole  Integrals 

The  residues  given  in  Table  1  may  also  be  derived  by  applying  a 
theorem  due  to  Maxwell  for  the  surface  integral  of  the  product  of  two 
harmonics  of  the  same  degree;  see  Hobson  (1931,  p.  157).  This 
alternative  procedure  has  the  advantage  that  the  singularity  is  located 
at  the  center  of  the  sphere.  It  will  be  illustrated  by  rederiving 
residue  ( 4) . 

Associated  with  the  multipole  potential  (1)  of  degree  -(q+1)  we  have 
the  regular  homogeneous  harmonic  of  degree  q, 

2q+1  ,  v 

Q  =  R  (40) 

q  P 


Then 


u  =  V< $  = 

P  P 


2q+1 


R 


[vaq  - 


(2q+1 )aq  p 


The  integral  in  (4)  then  becomes 


(41) 


/  n •Up'FdS  =  /  tn.Vfi  -  £  (2q+1)ft  IMS 

S  c  S 

■  -  I  Vds  *  -  "iiTJ  /  Ws  1421 

c  S  c  S 

by  (17)  and  (16a),  and  the  orthogonality  of  the  harmonics.  The  Maxwell 

formula  then  gives 


/  VqdS 

s 


4ttM  2q+2 
2q+1  C 


s  Tq 


(43) 


in  which,  in 
Substitution  of 


the  limit  as  c-x>,  'F  may  be  replaced  by  by 

q 

(43)  into  (42)  then  gives  the  residue  (4). 


(21  )  . 
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